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Abstract 

We determine the next-to-leading order dispersion laws for slow-moving quarks in hard-thermal- 
loop perturbation of high-temperature QCD where weak coupling is assumed. Real-time formal¬ 
ism is used. The next-to-leading order quark self-energy is written in terms of three and four 
HTL-dressed vertex functions. The hard thermal loops contributing to these vertex functions are 
calculated ab initio and expressed using the Feynman parametrization which allows the calculation 
of the solid-angle integrals involved. We use a prototype of the resulting integrals to indicate how 
finite results are obtained in the limit of vanishing regularizer. 
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I. INTRODUCTION 


The past ten years or so have witnessed an abundant activity that tries to understand 
the properties of the quark gluon plasma, the mechanism(s) of deconhnement and the char¬ 
acteristics of the transition from hadronic matter to quarks and gluons. Experimentally, 
collaborations at RHIC 111 and now ALICE at the LHC yf are among the main efforts 

n 

dedicated to this aim. On the other hand, lattice simulations [3[ as well as hydrodynamic 
modeling jd] help investigate the thermodynamic and transport properties of the plasma. 

From a perturbative QCD standpoint, calculations at high temperature use the so-called 
hard-thermal-loop (HTL) summation of Feynman diagrams j^. For example, one deter¬ 
mines the pressure and the quark number susceptibilities from the thermodynamic potential 


calculatec 


plasma j?!. 

HTL summation came about in order to overcome early problems encountered in the stan¬ 
dard loop-expansion of high-temperature QCD 8j. However, it makes the next-to-leading 
order (NLO) dispersion relations for slow-moving quasiparticles, quarks and gluons, difficult 
to calculate as they involve the use of the fully HTL-dressed propagators and vertices. The 
first quantity calculated in NLO fully-HTL-dressed perturbation is the non-moving gluon 


to two and three-loop order j^, or the electric and magnetic properties of the 


pte H. That waa followed by the calculation of the nou-.uovuig quark daurp- 
ing rate^ These calculations have been performed in the imaginary-time formalism of 


hnite-temperature quantum field theory 12[; they extract the damping rates from the imag¬ 
inary part of the fully-HTL-dressed one-loop order self-energies after analytic continuation 
to real energies is taken. In a line of works, we have used this formalism and looked into the 
infrared behavior of fully-HTL-dressed one-loomorder damping rates of slow-moving longi- 


and quasiparticles 


tudinal 111 and transverse gluons [l^, quarks [l^, fermions^ [l^ and photons [l7| in QED, 




181 in scalar QED^. 


In this logic, the natural step forward is to try to calculate the NLO energies of the 
quasiparticles. That would come from the real part of the fully-HTL-dressed one-loop order 


^ This quantity was later calculated in Uj using the real-time formalism. 

^ Note that the fermion damping rate at zero momentum in finite-temperature QED is independently 
calculated in 


llj using the real-time formalism. The same result is found. 


In scalar QED, we have also calculated the NLO energy of the quasiparticle. 
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self-energies. This is notoriously much harder than extracting the imaginary part. The hrst 
contribution in this direction is the determination of the pure-gluon plasma frequency Ug (0) 
at next-to-leading order in the long wavelength limit 19| . Imaginary-time formalism is used 
and the number iVj of quark flavors is set to zero from the outset. A gauge-invariant result 
is found: 

VK 


oog (0) = ^gT (l - + ...), Nf = 0. 


( 1 . 1 ) 


In this result, g is the strong coupling constant, T the temperature and Nc the number of 
colors. The next contribution came some time later [^, [^| and 21|, namely the determi¬ 
nation of the NLO fermion mass ujg (0) in high-T QCD (and QED). For quarks, the result 
found is |2l| : 

a.,(0) = ^('l + ^9 + ...j, N, = 3,N, = 2. (1.2) 

This calculation was performed in the real-time formalism of quantum held theory (for 
reviews on this formalism, see 22|b 

One should note in this respect that the NLO contributions to such quantities come for 
soft one-loop diagrams. Indeed, a general power-counting analysis performed in 23| using the 


real-time formalism shows that, except for the photon self-energy where two-loop diagrams 
with hard internal momenta do contribute, next-to-leading order contributions come from 
soft one-loop diagrams with HTL-dressed vertices and propagators. The work [h] shows 
that the usual power-counting in imaginary-time formalism overestimates a number of terms 
that are in effect subleading. 

The present work aims at determining the NLO dispersion relations, real part (energy) 
and imaginary part (damping rate), for slow-moving quasi-quarks in a quark-gluon plasma 
at high temperature with bare masses taken to zero. We use the closed-time-path formu¬ 
lation of the real-time formalism of hnite-temperature quantum held theory |2J, |2^ . The 
advantage is that we avoid the analytic continuation from discrete Matsubara frequencies to 
continuous real energies and all that comes with it which, in a sophisticated calculation like 
the determination of the dispersion relations, can make it difficult to extract the analytic 
behavior of the physical quantities. But as everything comes with a price, one disadvantage 
is that, as a result of the so-called doubling of degrees of freedom, each n-point function 
acquires a tensor structure with 2"' components to start with, which means a signihcant 
increase in the number of say one-loop diagrams involving three and four-point IPI vertex 
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functions. In addition, this calculation will not benefit from nice simplifications that arise 
when we set the quark momentum to zero, like the replacement of momentum contractions 
of HTL vertices with appropriate HTL self-energy differences via Ward identities 21|. 

This article is organized as follows. After this introduction, we dehne in section two the 
HTL-dressed quark and gluon propagators, as well as the quark energies and damping rates 
at next-to-leading order g'^T. These quantities are directly related to the NLO fully-HTL- 
dressed quark self-energy Sd). This quantity is calculated in section three. We give there an 
explicit expression of Sd) in terms of the three and four HTL-dressed vertex functions. These 
functions are derived ab initio as discrepancies between different results in the literature are 
found 2m, l27|. Then, in section four, we introduce a Feynman parametrization to help 


perform the solid-angle integrals present in the vertex hard thermal loops. 

Still, the subsequent integration task remains formidable. In section hve, we take a 
prototype and show how one can carry out with such integrals. The work is mainly numerical. 
We choose to avoid using the spectral decompositions of the HTL-dressed propagators and 
aim at getting a hnite result with the multi-integral as dehned. We indicate in this section 
how it is possible to obtain a stable behavior down to 10 “® in unit of the quark thermal 
mass nif. 

Brief concluding remarks populate section six. An appendix is dedicated to the derivation 
of the three and four-vertex hard thermal loops. 


II. THE NLO DISPERSION RELATIONS 

We consider QCD with colors and Nf flavors. The quark dispersion relations can been 
cast as: 

det (f-E (P)) = 0 . ( 2 . 1 ) 

Here, P = {po,i^ is the quark external soft four-momentum and E (P) is the quark self¬ 
energy, which can be decomposed into two components Ej. in the following manner: 

E (P) = 7+pE-m + 7-pS+(i’). (2.2) 

In this expression, = ( 7 ° T 7 -p) /2, with p = p/p and { 7 ^} the four Dirac matrices. 
Relation fl2.ip is equivalent to the following two dispersion relations: 

Po TP - (P) = 0. (2.3) 
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On shell, the (complex) quark energy pq = Q (p) can be decomposed in powers of the coupling 
constant g: 

n (p) = (p) + (p) + ... . (2.4) 

This follows a similar decomposition of S, namely: 

S(P) = 

Shtl (P) + E(i)(P) + ..., (2.5) 

where Shtl is the lowest-order contribution, formed by the hard thermal loops of order gT, 
and the NLO contribution, of order g^T. The contribution (p) = ijj± (p) is thus of 
lowest order pT, and (p) is the NLO contribution of order g‘^T . The dispersion relations 
fl2.3p can therefore be decomposed as: 

uj±{p) + (p) H-= ±p + Shtl± (f^± (p),p) + (0± (p),p) + ... . (2.6) 


Remembering that p ~ pT, We have: 

S2St"(p) = 


s±’ (‘i'±(p).p) 


(2.7) 


I - ShTL± (^,P)L=a;±(p) 

Here d^i stands for d/du. The real parts of o£^(p) are the NLO corrections (p) to the 
plasma quark energies, and the negatives of the imaginary parts are their damping rates 
7± (P)- 

The quantities uj± (p) are the solutions to the lowest-order dispersion relations in which 
only Shtl± {P) are retained in (12.3p . These latter are known: 


^HTL± 


{^,P) = 


m 


2 r 


P 


±1 + :^ It 


LJ 


In 


UJ +p 


( 2 . 8 ) 


2 \ p / 00 — p 

where ruf = ^/Cf^gT is the quark thermal mass to lowest order with Cp = {N^ — 1) /2W- 
The lowest-order quark energies u± (p) are real; they are displayed in Fig. [H Note how 
quickly the ultra-relativistic behavior sets in, at already p ~ 6(2)m/ for a;+(_)(p). This 
indicates that the soft region is effectively narrow. For soft p = p/mj < 1, they can be 
obtained in power series: 


( 1 1 16 
1 ± gP + T + O (p^) 


Also, from the dehnition of the HTL self-energies Fhtl, oae can rewrite: 


^±\p) = 


ul (p) 


2mj 


—(‘^±(p).p) 


(2.9) 


( 2 . 10 ) 
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FIG. 1: Lowest-order quark energies. Ultra-relativistic behavior sets in quickly. 


The HTL self-energies Shtl define also the HTL-dressed qnark propagator, which can 
also be decomposed into two components: 


A (P) = T+, A_ (P) + 7_,A+ (P) ; 

,2 r 1 


(P) = po±p- 


m 


f 


p 


Po+P 


=Fl -k —mi {p ± po) In 

2p ^ Po- p 


( 2 . 11 ) 


The HTL-dressed glnon propagator is also a qnantity we need. In the Landan gange^, it is 
given by the following relation: 

A„„ (K) = At (K) Pi + At (K) (2.12) 

in which P^^ are the usnal transverse and longitndinal projectors respectively: 


^ai/ ~ 9^lV 






= 

flip 


K„K, 


(2.13) 


^2 ^2 ’ Ml " ^2 ’ 

where, in the plasma rest-frame, K = (^k, kok^. The qnantities A(r ^ are the transverse and 
longitudinal gluon HTL-dressed propagators respectively, given by: 


A^^{K) = K^-m; 


' ^ K'^ /co , ko + k 


AZ" (K) = 




ko ko + k 

1--In --- 

2 k ko — k 


(2.14) 


^ The Landau gauge is part of a class of covariant gauges for which the soft one-loop order corrections to 
the lowest-order dispersion relations are independent of the gauge |19|. 
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In this expression, rUg = {Nc + Nf /2)/6gT is the gluon thermal mass to lowest order. 


III. THE NLO QUARK SELF-ENERGY 


There are two one-l oop HTL-dressed diagrams that contribute to the NLO quark self¬ 
energy displayed [28| in Fig. [2] and Fig. [3l The diagram in Fig. [2] writes as follows: 



FIG. 2: NLO HTL-summed quark self-energy The large dots indicate HTL-dressed vertex 

functions and propagators. All momenta are soft. 


sf’(C) =-i/Cp [ ^T>‘(P,Q)A(Q)T‘‘(Q,P)A,,(K), (3.1) 

J (Zir) 


(2 

with Q = P — K, and the diagram in Fig. [3] writes as: 


hi) 


^2 iP) = 


-ig^C^ 


d^K 


F^" (P,K) {K). 


(3.2) 


There are three summation structures: Lorentz (explicit), Dirac, and RTF. We introduce 


now the Keldysh indices (“r/a” basis) of 


finite-temperature real-time formalism 22 


he closed-time-path (CTP) formulation of the 


25|. The retarted (R), advanced (A), and sym¬ 


metric (S) propagators are given by the following definitions: 


(K) = A^bV (^) = ^b,f [ko + ze, k) ; 

A^^^ (K) = A”p (K) = /\b,f [ko - te, k) ; 

(K) ^ AIb (K) = Nb,f (ko) sign(fco) [A|,^ (K) - A^b (K)] , (3.3) 


where B stands for bosons and F for fermions, and Nb,f are related to the Bose-Einstein 
Fermi-Dirac distributions ub^f via the relations: 


Nb,f (x) = 1 ± 2nB,F (a:); nB,F (x) 


1 

ehl/'T =p 1 ‘ 


(3.4) 
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We then have for the two components of the following explicit expressions: 

+r^„ (P, Q) (Q) rr„ (Q, P) a;„ (P) 

+r^„ (P Q) A® (Q) r:„ (Q, P) A,t (P) 

+r;;„ (p q) a" (q) r:„ (q, p) a;j„ (p) 

+ r'‘„ (P,Q)A^(Q)T-^{Q,P)A^(K)], (3.5) 

and for the two components of the following expressions: 



( 2 ) 

FIG. 3: NLO HTL-summed quark self-energy . The large dots indicate HTL-dressed vertex 
functions and propagators. All momenta are soft. 


(C) = 


-ifcp r cPk 

~ J (it)* 


ti-7±p[rr„(CP)AS„(p) 


+ rc, (P P) A« (P) + rs;, (P p) a* (p)] , 


(3,6) 


Note that these expressions of s[2 in Eq. (13.51) and in Eq. (13.61) have been written in 
ll| using a different notation. 

29| . However, 


The HTL-dressed vertex functions T are derived in the literature 
there are discrepancies in these results 


26 


26 


m 


27j, which led us to rederive all three and four- 


point HTL-dressed vertex functions ab initio in the CTP formalism. We recover the results 
of [27j; these are presented in appendix A. For our needs, we have: 


r;;„(PQ) = r + i"-(P,Q)-, 
n„(pQ) = 7 '‘ + /(;-(pq); 

K.AP,Q) = K,JP,Q) = 0, (3.7) 
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for the two-quarks-one-gluon vertices, and: 

(P, K) = /r_ (P K ); TZ, (P A') = r;;, (p A) = 0, (3.8) 

for the two-quarks-two-gluons vertices. In these relations, the quantities J’s are hard thermal 
loops given by these solid-angle integrals: 


y.,yA,(3) = m} 






An {PS + irj^e) {QS + iri2e) ’ 

(PK) = f 

^ ^ f J An {PS + t7]^e) {PS + irj^e) 

1 1 


X 


_{P + K)S + i7]^e ' {P-K)S + i7]2e\ ’ 
where S = (1, s) and the indices r]i and 772 take the values -|- or —. Using these results, we 
can rewrite and ^22 above as simply: 


d±(A) = 


-ig^CF f dfK 

2 J 


[(7'* + F__ (P,Q)) (Q) (Y + r_ (0, p)) a;;„ (k) 


E«(F) = 


(Y + A- (P «)) (0) (Y + II- (Q, P)) (A); 

d^K 

trj^P’i{P,K)Al{K). 




(3.10) 

n4 i±p" - 1 I \^^ / ■ (3.11) 

^ J {2n) 

We now work out the contractions over the Dirac and Lorentz indices. Starting with 
the contributions that do not include a hard thermal loop write generically as: 

(P A) = tr (o-.p-^Pp") Al (A) A!',, (Q) 

= -2 ( 1 -p.-M.-b Aj(A') AY,. (<3) 

2/co/c 

X i^p,.k + qe.k^ +k'^{l+ pe.qe) {K) {Q) . (3.12) 

In this expression, Ap {K) = Ap {K) /. Also, p^ = with = ± and similarly for g^, 
with summation understood over e^. The superscripts X and Y indicate the RTF indices 
(R, S, and A). The contributions that involve one hard thermal loop vertex function writes: 


kill 


p,.q, - 2p,.kq,.k 


tpXY 

X 


(P A) = tr (7..4y.7..T‘') A;;„ (A) A':,, (Q) = m} 


dQg 


1 — pe-qe ~ Pe-S — (Je-S + pe-k k.S -|- q^.k k.S — k.S^ + pF^.qek.S^ + 2p^.S (Je-S 
-p^.kqe.sk.s - qe.kp^.sk.s^ Ap {K) AI^^ (Q) - {k^ {1 + p^.q^ - Pe-S - qe-s) 
—kok (pe-k + qe-k + 2p^.q^k.s — q^-kp^-s — p^-kq^.s — q^.s k.s — p^.s k.s 


An {PS + ipi£) {QS + ip 2 Q 

(3.13) 


Pk^ipe-k k.s -|- pe-k k.s — k.s -|- p^.pek.S — pe.kpF.S k.s — pe.kpe.S k.s ) ) Ap {K) A_^ {Q) 
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Here 77 ^ and 772 ^^ke the values ±. Due to the symmetry in The other contribution 
with one hard-thermal-loop vertex function is equal to the one above when changing 77 ^ into 
772 , namely: 


tr («) = t"- (^0 Ay, (Q). (3.14) 


The contribution involving two hard-thermal-loop vertex functions is longer. It generically 
writes: 


'^p-,VlV2\VlV2 


(P,K) ^ tr (7,yy„7,.7,y;) A^ (A') Ay, (Q) 

4 r dQ.. 1 


= m 


^ J in {PS + ir]^e) {QS + irj^e) 


X 


dVL.i 


- -[(—s.s' — p^.q^s-s' + pe-s s.s' + q^.s s.s' + p^.s' s.s' 

_L rtrn' C- \ 


An {PS' + 7772^) {QS' + 777^y) 

+qe-s' s.s' — pe-s q^.s' s.s' — pe.spe.s' s.s' — s.s'"^ + p^.qeS.s''^ + k.s k.s' + p^.q^^.s k.s' 

—pe.s k.s k.s' — q^.s k.s k.s' — p^.s'k.s k.s' — q^.s'k.s k.s' + p^.s q^.s'k.s k.s' 

+q^.sp^.s'k.s k.s' + s.s'k.s k.s' — p^.q^s.s'k.s k.s''^ {K) A^^^ {Q) 

- (1 -h pe-F - Pe-s - qe-s - pe-s' - qe-s' + Pe-S q^-s' + qe.sp^.s' F s.s' - p^.q^s.s') 

—kok {k.s + k.s' + pe.qe k.s + pe-F k.s' — p^.s k.s — pe.s k.s' — qe.s k.s — qe.s k.s' 

—pe.s' k.s — pe.s' k.s' — (Je-S' k.S — (Je-S' k.s' -|- pe.S (Je.s' k.S -|- pe.S (Je-S' k.s' 

Pqe.spe.s' k.s + qe.s Pe.s' k.s' + k.S S.s' + k.s' s.s' — pe.qek.S s.s' — pe.qek.s' s.s'^ (3.15) 

-y/Cg {^k..S k.s' + pe.qek.s k.s' — pe.s k.s k.s' — pe.s k.s k.s' — Pe.s' k.s k.s' — Pe.s' k.s k.s' 

+Pe-S Pe.s' k.s k.s' + pe.Spe.s' k.S k.s' -f k.S k.s' S.s' + pe.pek.S k.s' S.s'j A^ {K) A^^^ {Q) . 

The integrand in is faster to write: 

1. / \ . Y . o /" dQ. 1 


^p',VlV2 


{P,K) ^ 2^r(7,/,7,^)Aj,(if)=7nJ 


An [PS* -1- 777;^y] [PS + ip 2 ^] 


X 


(P -|- P) S' -|- ipiE {P — K) S + 7772 ^ 

X — Pe.s — k.s^ + Pe.s k.s^^ Ay {K) (3.16) 

k^ (1 — pe-s) — 2kQk {k.s — pe.s k.s) + k^ pe.s k.s'^') Ay {K). 


From these expressions, we can write the NLO HTL-dresses quark self-energy in a compact 
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form: 


E™ (P) = 


-ig^Cp f d^K 


[F|5 (P, K) + F^l (P, K) + 2F^^_ (P, K) + Fg__ (P, K) 


2 J (27r) 

+Fg-+ (P, K) + p|^__.__ (P, i^) + (^, K) + Gl... (P, 7^)] . (3.17) 


IV. HTL VERTEX FUNCTIONS AND FEYNMAN PARAMETRIZATION 


The next step is to find a way to evalnate the solid-angle integrals involved in the hard- 
thermal-loop vertex fnnctions. The way we do this is to rewrite these integrals nsing the 
Feynman parametrization. From Eqs. fl3.13p . fl3.15p . and fl3.16p above, we see that we have 
two kinds of solid-angle integrals to deal with, namely: 

dVL, S>^S^ 


= 

CZ (-P. K) = 


dvr [PS* ir]^e] [QS* iri 2 £] ’ 
d^s 

dvr [PS' -|- iriis] [PS* -|- iri 2 ^] 

1 1 


X 


+ 


_{P + K) S + {P-K)S + ir]2£ 

The simplest of all these integrals is the integral: 


(4.1) 


(4.2) 


dvr {PS + irjis) {QS + ir] 2 £) 

Remember that S is the d-vector (1, s) and the integration is over the solid angle of s. Let 
ns pnt aside the ie prescription for the moment. Using the Feynman parametrization: 


AB 


= / du- 


1 


we write: 


7-00 

'^VlV2 


{P,Q)= / du 

Jo 


lo {Au + B {1 - u))'^' 

fdQs 1 _ 

dvr [{P-Ku)Sf ~ Jo {P-Kuf 


du 


where K = P — Q. In this case, the integral over u can be done formally to get: 




(1 - Ml) U2 

2\/A (1 - M 2 ) Ml’ 


In 


(4.3) 


(4.4) 


(4.5) 


with the notation ui ^2 = [PK ± y/Aj and A = PK^—P^K^. Noticing that PS+irhe = 
Po + 'ipp ~ P-P the reintrodnction of the ie^s in the final resnlt (ld.5p is a matter of shifting 
Po ^ Po + ipp and go —^ <?o + This also applies to the two next integrals. 
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The next HTL vertex function to consider is the solid-angle integral: 


j 4^ (p5 + ^Qs + ^r^^e) • 

Using a Feynman parametrization and the notation R = P — {P — Q)u in which po and Qq 
are redehned with the corresponding fe’s, we obtain the result: 




Li (ro — xr) 

0 1 ro r r* 


"o — r^ 2r ro — r J r^ 


Little useful comes from pushing further the integration over u as the hnal result will not 
have a reasonably simple formal expression; it is better to leave this expression as it is and 
let the integration over u be done numerically. However, the Feynman parametrization is 
useful as it reduces the number of integrations to be performed from two (solid angle) to 
one (over u) or zero when this latter is explicitly feasible. 

The third HTL-vertex solid-angle integral is: 


J'= (P, Q) = 


An {PS -|- irj^e) {QS + 


This solid-angle integral can be performed using the Feynman parametrization and, still 
with R = P — {P — Q) u, we have: 

Jo 

^^1^2 2 1^ O ^ ) ’ 

y Ir Vq — r J 


^VlV2 


_ ^2 .y ^2 


ro To+ r 

1-In- 

2r Tq — r 


The integration over u is left to be performed numerically. 

The two-gluon-two-fermion HTL vertex function uses the Feynman integral: 
1 

■ ^ — 2! / dui I du2 -Q- 

ABC Jo Jo [^C-B)uiU2 + {B-A)ui + Af 

Now calling Pi 2 = (po + rewrite: 


(4.10) 


1,7“ (P, K) = 

JiJZ K) = 


Ji;z (P, K) + JZZ (P, -K ); 

r (Pis s^^s•'s^ 

J An P 1 SP 2 S {Pi + K)S' 


(4.11) 
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Calling C = (Pi + K) S, B = PiS, A = P 2 S, and the four-vector T = U 1 U 2 K+U 1 (Pi — P2) + 
P2, with to = uiU 2 ko + iui (77^ — 773) e + po + t = uiU 2 k + p, we have: 

dQs S'^S'^S'^ 


J,7“(P,P) = 2/ dm 7/1 


dM2 


(to-t.))' 


(4.12) 


Since the quantity is fully symmetric in the Lorentz indices, we only need to hnd the 
000, OOi, Oij, and ijk components. We have: 

dhl, 1 


= 2 / dmm 


duo 


= 2 


dui Ui / du2 


(to - t(s)^ 

to 


(4.13) 


.70 .70 {to-t‘^)‘^ 

Note that in this case, the Feynman parametrization does not reduce the number of integra¬ 
tions to perform. Yet, it still has an advantage over the original integration over the solid 
angle of the original unit vector s as, cast this way, the subsequent integration over the 
azimuthal angle of k will now be trivial. The next solid-angle integral is: 


jJ°;(P,P) = 2 / dm 771 


dU2 


dfl. 


{to - t.sy 


= 2 


dui ui / du 


0 - (‘Y 


(4.14) 


The next solid-angle integral involves a symmetric rank-2 tensor structure: 


J^%^{P,K) = 2 / dm77i 


duo 


dhlo 




S S 


= 2 


'0 


477 (to - t.s) 

dui Ui I du 2 ; 

to +t 


'0 


to 


2t2(t2-t2) 4t3“^to-t’ 


iln 


Dr^ 


to ( 5 t^ - 3 to) ^Y^to + t 


(4.15) 


2«n«§-e)^ ' 4i3‘“«„-i’ 
and the last solid-angle integral involves a completely symmetric rank-3 tensor structure: 


jyXiPi^) = 2 / dUYUY 


duo 


dklg 


>0 


'0 


= 2 


47r {to-Uy 

dui 771 [ dU 2 {fSt^ + ■ 


-^1)11)2 2t3 ’ ^ 4" 


F = 

^ V 1 V 2 


''0 


3fo , ^0 + ^ 

— In 


(tg-P) 2 t to-tj' 
5 


2 + 


P 

''0 


3p , to + t 

-In 


{tl - pf 2 P V (to - ty 2 t to-t 


(4.16) 
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V. INTEGRATION: A PROTOTYPE 


The next task is of course to evaluate all the integrals in Eq. fl3.17p . one after the other, 
and add them together. This is not an easy matter. Here we show how one can carry out 
with such expressions. We work in units of the quark thermal mass, i.e., we take ruf = 1. 
This means that the ratio: 


m 

m 


16W(W + iV;/2) 

6(N;-1) 


= 3 + 


N, 


(6.1) 


with Nc = 3 will replace in the gluon propagators. This ratio is equal to 4 for Nf = 2, a. 
value we take in this section. 

It is best to work out a prototype, namely®: 


/ (t, e) = I d^KA% {K) (g) 


dVL. 


dvr {PS — is) {QS — is) 

f'+OO ^+00 nl 

dk dx dko / du M (t, k, k^, x, u, e) 

) J-l J-oo Jo 


(5.2) 


The notation is explained as we proceed. As before, P = {po,p) is the four-momentum of 
the external soft ”on-shell” quark, K = {ko, k) the four-momentum of the soft gluon in the 
loop, and Q = P — K] see Fig. |2J The quantity x is cosine the angle between k and p - the 
integration over the azimuthal angle of k is done trivially. The integration variable u is from 
the Feynman parametrization when the solid-angle integral is performed, see Eq. fl4.4l) . The 
quantity t is the ratio p/po, to be introduced more naturally later in this discussion. The 
integrand M in this example is therefore given by the following relation: 

As (/co) sign (fco) r.R 


M {t, k, ko, X, u, e) = 


{k, ko, e) — Ay {k, ko, {q, go, e). (5.3) 


-|- 2Bu + A 

The quantity Nb {ko) is related to the Bose-Einstein distribution, see Eq. fl3.4p . and 
A^{k,ko,e) is the retarded transverse gluon propagator, which can be written as: 


A, 


R(-l) 


{k,ko,e) = ^ - kl) 


' \ ((^0 -k)'^ + e 

3P ^\{ko + k)‘^ + e^ 


2\ H 


(5.4) 


+* 


2ko 


(/cg — /c^) I arctan 




ko — k 


— arctan 


ko P k 


- esign(/co) 


® Prefactors are not important in this section. They are not displayed. The quantities Re/(t, e) and 
Im/(t,e) are plotted in arbitrary units. 
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The quantity A^' {q, qo, s) is the retarted quark propagator, given by: 




{q,qo,£) 


1 

- + q-qo- 
q 


qo-q 

4g2 


In 


f{qo±ql+f\ 

\{qo - qy + £^) 



arctan 


qo + q 


— arctan 



—i 



+ 


qo-q 

2g2 


arctan 


{qo±ql±^\ 

{qo-qy + e^J 

( -) — arctan 

\qo + qJ 



(5.5) 


The momentum q = p — k and the energy qo = po — ko where po is the quark energy. The 
quantities F, B and A in the denominator of M are regularized relativistic relations: 


A = (po- ief - 
B = k.p — {ko — 2ie) {po — is ); 

F = {ko- 2ief - k^ (5.6) 


As mentioned above, the variable t is the ratio p/po, in terms of which the quark momentum 
writes: 

p (*) /’"/ = 

Remember we take the quark thermal mass mf = 1. The above relation comes from the 
lowest-order quark dispersion relation Al^(po,p) = 0 where Al^ is given in Eq. fl2.1ip . The 
behavior of p{t) and po{t) = oj-{p{t)) is shown in Fig. 01 Slow-moving quarks would have 
p < 1 (in units of mf), which translates into t < 0.64. The limit t —)■ 1 is the fast-moving 
region in which both p{t) and po(t) are large. 

In this example, the integral over u can be done analytically, which reduces the number 
of integrations in Eq. fl5.2p to three. Indeed, one obtains: 

^+oo ^+1 r+cxD 

I{t,e) = dk dx dko k'^N b {ko) sign {ko) 

Jo J-l J-oo 

X [A^ {k, ko, e) - A^ {k, ko, -e)] {q, qo, e) V {t, k, ko, x, e), (5.8) 
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FIG. 4: The quark momentum p (solid) and energy pQ (dashed) as functions of the ratio t = p/po- 
The quark mass ruf is taken equal to one. The fast-moving limit is t —?• 1. Slow-moving quarks 
would be for t < 0.64. 


with the function V given by the following expression, see Eq. 


V (t,k,ko,x,e) = / du ---- 

^ ^ Jo Fu^ + 2Bu + A 

2\/A I 2 


((ri - 1)^ (r2^ + i2^) 


((r2-l) +i2^) (rf + ii^) 


+i 


arctan I — 1 — arctan ^ ^ 
*1 


arctan I — ) arctan 
A2 


r2 - 1 
12 


(6.9) 


The quantities are = Re Sj and f, = Im Si with: 


Si = 


-R + ^/A 

F ’ 


S2 = 




A = _ af. 


(5.10) 


There are difficulties with the integral I {t,e). One is that the integrand features sudden 
jumps, those coming from the arctans in the propagators, the Bose-Einstein distribution at 
ko = 0, and eventually those coming from V. Such jumps make any integration method 
either too long to be useful, or give unstable results. This instability is more intense when 
trying to take e smaller and smaller. So, in order to see more closely what is at stake, we 
have partitioned the integration region in the {k, ko) plane into domains bounded by the 
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following lines of sudden jumps: 


/co = 0; /co = ±fc; fco = Po ± \/— 2pkx; 

k = k,^ = ll^J^-llnfl±^y 

2po — xp 2t 1 — xt]J 1 — t 2 \1 — t J 


(5.11) 


The last (vertical) line is simply the line k = pa—q. Using the change of variables 9 = arctan k 
and (f) = arctan fco, these domains are plotted in Fig. |5l If one integrates M in each of the 



FIG. 5: Boundaries in (fe, /cq) plane at which the integrand M has sharp jumps (see text for details). 
Here t = 0.45 and x = 0.8. 


domains of Fig. |5] separately and sum up the results, one hnds that / {t,e) becomes stable 
and reliable in the limit £ ^ 0. See Fig. |6] (the units of I{t,e) there are arbitrary): both 
the real and imaginary parts of I{t) behave smoothly, and these two plots are obtained with 
£ = 3.2 X 10“® (in units of m/). 

The dependence in £ is of course an issue to explore. Fig. [7] shows how for example 
I{t = 0.32, e) behaves as a function of m = — logj^gS. The behavior is stable down until 
£ ~ 10“® (in units of ruf), below which the numerics loose reliability. This is way beyond 
any precision we can hope for, and we notice that I{t, e) converges smoothly to a hnite value 
already satisfactorily reached at £ ~ 10“®. The smooth convergence to hnite values is also 
obtained for other values of t. 
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FIG. 6: The real and imaginary parts of I{t) for soft p (meaning t < 0.64). Both parts have a 
stable behavior. Here, e = 3.2 x 10“®. The units of I{t) are arbitrary. 




FIG. 7: The e-dependence of I{t). We see stability and convergence up until e ~ 10 Here 
t = p/pq is take equal to 0.32. 

VI. OUTLOOK 

The present work aims at calculating the energies and damping rates of slow-moving 
quarks at next-to-leading order in the context of massless QCD at high temperature using 
the real-time formulation of the fully-dressed hard-thermal-loop perturbative expansion. 
These quantities are extracted from the poles of the quark propagators. At lowest order, the 
energies u± {p) are real, see Eq. fl2.9p . The next-to-leading order contributions necessitate 
the determination of the next-to-leading order quark self-energies see Eq. fl2.7p . 

In this work, we have given the analytic expression of in terms of loop-four- 
momentum integrals involving fully-HTL-dressed quark and gluon propagators and vertex 


18 











functions, see Eqs. fl3.5p and ( I3.6p . These expressions were already written in [ll|. The 
HTL vertex functions themselves are derived and written as solid-angle integrals. We have 
rewritten these latter using the Feynman parametrization in order to perform these solid- 
angle integrals. 

The next step is to perform the integrals. This is done numerically. The usual approach 
is to use the spectral decomposition of the dressed propagators. But we refrain from doing 
so and will try to tackle these integrals head on. This is known to be difficult. One main 
difficulty is the jumps the integrands suffer from at the singular lines of the propagators, 
more prominently as the regularizer £ —)■ 0. We have used a prototype integral to indicate 
how such difficulties might be overcome. By splitting the integration region into appropriate 
domains, we can obtain an estimate of the prototype integral with a robust behavior down 
to e ~ 10 “® in unit of the quark thermal mass rrif. 

Of course there are many other terms to handle, more involved than the prototype. Will 
there be stability for each? Can we add them all? This is currently under investigation. 


Appendix A: Hard Thermal Loop Dressed Vertex Functions 


Discrepancies between different results in the literature 2g, l27| regarding the derivation of 


the three and four-point hard-thermal-loop vertex functions in the CTP formalism imposed 

n 

on us an ab initio recalculation of these quantities. We recover the results of 1271. This 

□ 

appendix summarizes our steps, with a notation close to [27J . 


1. Three-Point Vertex Functions 

The quark-gluon vertex functions are the sum of the bare vertices and the corresponding 
hard thermal loops: 

(Ai) 

In the {12} basis of the CTP formalism, the bare vertex 7 ^ is given by the relations: 


(— 1 )* ^ 7 ^ when i = j = k 
0 otherwise 


'yijk 


(A 2 ) 
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The indices i, j, and k take the values 1 and 2. The hard-thermal-loop contributions (5r^ to 
these vertices are obtained by calculating the following one-loop diagrams: 


ir;;, (p, q,r) = - 2 iy/v, / (y. (A 3 ) 

J 127rj J 127rj 

in which the loop momentum K is hard, in front of which external momenta are neglected. 
In this expression, 

A, {K) D,k {K - Q) Dk^ {K + P)- 

V:jk = A, (K) A, {K - Q) D,, (K + P). (A4) 


The functions Dij [K) are the four {12} components of the bare bosonic propagator: 

Dn (K) Du (a\ 

D21 (K) D22 (K)) 

( l+riBikp) _ nB(ko) e{-kp)+nB{kp) _ e{-kp)+nB(ko) 

K'^-\-ie K'^—ie K^-\-i£ K^—is 

0{ko)-\-nB{ko) 0{ko)-\-nB{ko) nBjkp) l^nBjko) 

K'^-\-is K‘^—i£ K‘^-\-i£ K‘^—i£ 


The quantity Aj {K) is Dij (K) with ub (/cq) replaced by —np {ko). We write these quantities 
in terms of advanced (A), retarded {R) and symmetric {F) propagators: 




where: 


Ai = ^{F + A + R)- Du = ^{F + A-R)- 
D 21 = ^{F-A + R)- D22 = ^{F-A-R), 


R{K) 

A{K) 

F{K) 


Ojko) 9{-ko) _ 

K"^ + ie K"^ — ie ’ 
0{-ko) 9{ko) 

K"^ + ie K"^ — ie ’ 

( 1 ± 2 «b,.(*^o))(^ 


le 




(A6) 


(A7) 


The vertex functions in the {RA} basis are linearly related to the same functions dehned 
in the {12} basis. For example, the three-point vertex function Tbaa is given by the linear 
relation: 

^RAA = Fill -|- F 112 -|- F 121 -|- F 122 . (A8) 
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This relation applies to the kernels V and V. Using the expressions flA4j) of the functions 
V and U', and the relations flA6D . we arrive at the following expression: 


Vraa — + R1R2R3 + -U1A2A3 + RiF 2 A^ + RiR 2 F^), 


(A9) 


where, for short, 1, 2 and 3 denote the arguments K, K — Q and K + P respectively. Next 
is to integrate over the internal momentum K. This is performed in two steps: hrst over ko, 
done using the residue theorem in the fco-complex plane. In this regard, the terms A 1 A 2 A 3 
and R 1 R 2 R 3 give zero contribution each. The three others give the following contributions: 


FAA 

RFA 

RRF 


+i 

Svr^ 

—i 

Svr^ 

—i 

Stt^ 



kdkNB^pik) 



kdkNB^pik) 



kdkNB,p{k) 


dQs 

dvr {PS — is) {QS + is) ’ 
dVL, 

dTT {QS + %e){{P + Q)S-iey 

dSt, _ ^ _ 

dTT {PS — is) {{P + Q) S — is) 


(AlO) 


Here, the quantities NB,p{k) are defined in Eq. fl3.4p and S = (1, s) is a time-like unit four- 
vector in which s is nothing but k/k. Remember that every external momentum is neglected 
in front of K. When summing the contributions together, some will add up to a product of 
propagator denominators multiplied by ie, so are dropped, and some will stand. In the end, 
we hnd: 


2 ^- dn 


dTpAA = f / kdk {hb {k) + Up {k)) 




27r2 


0 j dvr {PS — ie) {QS -|- ie) 


Now the integration over k can be done. Using the known results: 

/•H-oo ^ 

d.kkun (k) = - _ , , _ 

Q .In ^ ' 12 


jr.2 r+oo 2 

dkkuB {k) = —T^; / dkknp {k) = —T^, 


Jo u Jo 

we hnally hnd the expression for this hard thermal loop: 

dQ Sf^^ 


dTpAA {P,Q,R) = -mj 


4p {PS - ie) {QS + ie) ’ 


(All) 


(A12) 


(A13) 
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with mj = Cp/SgT. The other three-point HTL vertex functions are obtained following 
a similar procedure. The relations between vertices in the {RA} and {12} bases are: 


^ ARA = 

= Fill 

+ 

Fi 12 

+ 

F2II 

+ 

F212 

^AAR = 

= Till 

+ 

ri 2 i 

+ 

r 2 ii 

+ 

r 22 i 

Fi^RA = 

= Fill 

+ 

Fi 12 

+ 

F221 

+ 

F222 

^ RAR = 

= Fill 

+ 

Fi 21 

+ 

F212 

+ 

F222 

F^iJi? = 

= Fill 

+ 

F2II 

+ 

F122 

+ 

F222 

^RRR = 

= Fill 

+ 

F122 

+ 

F212 

+ 

F221 


(A14) 


with TAAA = 0 identically. For each of these vertex functions, one works out steps similar 
to the ones for T/j^a and one hnds: 


6 Tara{P,Q,R) 

STaar {P,Q,R) 

6 Trrr{P,Q,R) 

6rRRA{P,Q,R) 


2 f dQs 

^ J dvr (PS -|- ie) (QS — is) ’ 

^ r dVts 

J An (PS ie) (QS ie) ’ 

2 f dn, S>^S 

^ J An (PS — ie) (QS — ie) ’ 

(P, Q, R) = 6Tarr (P, Q, r) = 0. 


(A15) 


2. Four-Point Vertex Functions 

The two-gluon-two-quark vertex functions are all hard thermal loops, no bare terms, 
given in the {12} basis by the relation®: 

d^K 




i+j+k+l 


(2n) 




X [CrD,, (K) Djk (K - Q) Dki (K + P + U) Du (K + P) 
-N,Di, (K) Djk (K - Q) Dui (K + P + U) Du (K + P) 


+ -N,Du(K - P)D,i(K -P- U)Du,(K + R)D,k(K) 


(A16) 


Consider for example the {i?A}-component 5VraaA) given by: 


dT raaa — ^Fiiii -i- 5riii2 -i- (5rii2i -i- (5ri2ii -i- 5rii22 + 5ri2i2 + 5ri22i + 5ri222- (ai7) 


A typo in (3.22) of [27| is corrected. 
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Consider the term multiplying (7^ in Eq. flA16p and call it Using the decompositions 

in Eq. ( IA 6 p . we obtain, in a similar symbolic notation as in Eq. flAQp . the relation: 


^^\iAAA — 2-^1^2^3^4 + + 2-^l-^2-^3^4 + ^-^l-^2-R3-^4- (^18) 

The subscripts 1,2,3, and 4 stand for the momenta iC, K — Q, K + P + U, and K + P 
respectively. The two contributions IA 1 A 2 A 3 A 4 and li?ii? 2 -R 3 -R 4 to yield zero each 

in the fco-complex-plane integration and are not displayed explicitly in Eq. flAlSp . The 
different terms in Eq. flAlSp contribute as follows: 

/■+°° kdk f dQ, iNB{k)Sf^S'^^ 


FAAA = 
RFAA = 
RRFA = 
RRRF = 


'0 


'0 


'0 


'0 


IGtt^ 
kdk 
IGtt^ 
kdk 
IGtt^ 
kdk 
IGtt^ 


dvr ((2P + U)SF is) ((P - R) S + is) {PS + ie) ’ 
dVLs -iNB{k)K^^K^$ 

dvr ((2P + U) S + is) {{R + P + U) S - is) ((P + U)S-iey 

dfl. 


iNF{k)S^^S'^^ 


dvr ((P -P)S-ie){{R + P + U)S- ie) {RS + ie) ’ 


dfl. 


iNF{k)S^^S^$ 


dvr {PS + ie) {{P + U)S- ie) {RS + ie)' 


Adding these terms together yield the following result: 


5rl 


RAAA 


— kdk {nB{k)+nF{k))S>^S‘'^ 

io J 4:71 {PS - ie) {QS + ie) {{P + U) S - ie) 


(A19) 


(A 20 ) 


The integrations over k can now be done using Eq. flA12l ). The terms multiplying Nc and 
Nc/2 in Eq. (IA16P are worked out in a similar way; they cancel each other. We therefore 
have the hard-thermal-loop four-vertex function: 

dQ, Sf^S’^^ 


STraaa{P,Q,R,U) = rrif 


dvr {PS — ie) {QS + ie) 

1 1 


^ _((P + P) ^ - ie) {{P + R)S- te)_ ■ 

The other {RA} four-vertex hard thermal loops are worked out in a similar way. With 
d^AAAA = 0 and using the linear relationships: 

d^ARAA = <^riiii -i- (5riii2 -i- 5rii2i -i- (5r2iii -i- 5rii22 + 5r2ii2 + (5r2i2i + 5r2i22; 

dT aara = ^riiii -i- 5riii2 -i- ^ri2ii -i- 5r2iii -i- dTi2i2 + ^r2ii2 + 5r22ii + 5r22i2; 

dT aaar = ^riiii -i- 5rii2i -i- 5ri2ii -i- 5r2iii -i- (5ri22i + ^r2i2i + 5r22ii + ^r222i; 

d^RRAA = <^riiii -i- (5riii2 -i- ^rii2i -i- 5rii22 + ^r22ii + ^r22i2 + <^r222i + ^ 1 ^ 2222 ; 

d^BARA = f^riiii -|- 5riii2 -|- 5ri211 - 1 - 5ri212 + ^r2121 + ^r2122 + ^r2221 + ^r2222; 

d^RAAR = <^riiii -|- (5rii2l -|- 5ri2ll - 1 - 5ri221 + (5r2ll2 + 5r2l22 + <^ 1^2212 + <5r22225(A22) 
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one obtains the following results: 


6Taraa{P,Q,R,U) = rrij 






X 


Att {PS + ie) {QS — ie) 

1 1 


_{{P + U)S + ie) {{P + R)S + ie) 

5Taara{P,Q,R,U) =mj • wnQ^ M 

'' J An {PS + le) {QS + le) 

1 1 


X 


STaaar{P,Q,R,U) = ruf 


_{{P + U)S + ie) {{P + R)S- i£) J ’ 

dSls 


An {PS + ie) {QS + ie) 

1 1 


X 


{{P + U)S-ie) {{P + R)S + ie)\' 

5Trraa {P, Q, R, U) = 6Trara {P, Q, R, U) = 5TnAAR {P, Q, R, U) = 0. (A23) 

The remaining eight components can be either calculated directly or obtained from those 
above using the KMS conditions. 


3. Change of Notation 


Finally, as we mentioned early in this appendix, the notation we use here for the vertex 
functions is close to that used in . However, the notation we use in the main text is close 


to the one used in 


21| . They are related in the following manner: 


= TP^,,{Pi,-P 2 ); 

^ 2 , Ps, Pa) = (A, P4, Ps, -P 2 ), (A24) 

with the understanding Ij = R{A) ■(-)■ ij = a(r). Thus, for the three-vertex functions: 
rj„ (P, Q) = (p -Q, P) = 7 " + p_ (P Q); 

r?„ (p Q) = {P,-Q,R) = r + R-(P,Q)-, 

rt,(PQ) = (P-Q, P) = 0 ; 

rf„ (p Q) = (P -(?, P) = 7 " + (P Q) ; 

r". (P Q) = (p -Q, p) = 0 , (A 26 ) 


with the following dehnition: 

/,yyp« = mj 


dn, si‘$ 

An {PS + iri^e) {QS + ir^^e) 


(A26) 
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For the four-vertex functions we need in the text, we have: 


TZ, (P, K) = TZ„ (P K,-K, P) = TZaaa iP-P -Q, Q) = 1-- (P K ); 

rc, (P K) = rz, iP K, -K, P) = tZaah (P -P -0. Q) = 

r^,“„ (P, K) = (P K,-K, P) = ^RARA (P-P -Q, Q) = 0. (A27) 


with the dehnition: 


(PK) = f 

J 4:71 [P5 -|- irjie] [PS* + f?72e] 

1 1 


X 


_{P + K) S + i7]^e {P - K) S + ir]2e_ ' 

Note that both and are already used in the main text, see Eq. (I3.9p . 
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